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Abstract 

Processes of heavy quark production at TEVATRON energies are considered using 
the semihard (/ct factorization) QCD approach with emphasis of the BFKL dynamics 
of gluon distributions. We investigate the dependence of the pT distribution of heavy 
quark production (presented in the form of integrated cross-sections) on different forms 
of the unintegrated gluon distribution. The theoretical results are compared with recent 
DO and CDF experimental data on beauty production. 
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1 Introduction 


It is very well known that from heavy (c— and b—) quark production processes one can obtain 
unique information on gluon structure function of the proton and the three gluon QCD vertex 
because of the dominance of the gluon-gluon fusion subprocess in the framework of QCD. 
These issues are important for physics at future colliders (such as LHC): many processes at 
these colliders will be determined by small x gluon distributions. 

Recently DO and CDF Collaborations have reported B.0 new experimental data on the 
cross section of inelastic open beauty production at the TEVATRON energies. A comparision 
of these results with NLO pQCD calculations shows that the calculations underestimate the 
cross section by a factor of ~ 2. Therefore, it would be certainly reasonable to try a different 
way. 

In the energy range of modern colliders (HERA, TEVATRON and LHC) the interaction 
dynamics is governed by the properties of parton distributions in the small x region. This 
domain is characterized by the double inequality s S> /i^ — s S> A^, which shows that the 
typical parton interaction scale /i (mass rric or px of heavy quark) is much higher than the 
QCD parameter A, but is much lower than the total center of mass energy y/s. The situation 
is therefore classihed as “semihard”: the processes occur in small x region (x ~ M'^/s -C 1) 
and the cross sections of heavy quark production processes are determined by the behavior 
of gluon distributions in this region. 

It is known also that in the small x region the factorization of QCD subprocess cross sec¬ 
tion and parton structure functions in a hadron assumed by the standard parton model is bro¬ 
ken 0-0. The resummation 0, H, ||, 0 of the terms [ln(/i^/A^) [ln(/i^/A^) ln(l/x) 

and [ln(l/x)as]” in the semihard (/cT’-factorization) approach (SHA) of QCD results in the 
unintegrated gluon distributions <h(x, q^, p), which determine the probability to hnd a gluon 
carrying the longitudinal momentum fraction x and transverse momentum qt at the probing 
scale p?. They obey the BFKL equation 0 and reduce to the conventional gluon density 
once the qr dependence is integrated out (/i^ = Q^): 

/ ^x,ql,Ql) dq^ = xG{x,Q^). (1) 

Jo 

To calculate the cross section of a physical process, the unintegrated functions have 
to be convoluted with off-mass shell matrix elements corresponding to the relevant partonic 
subprocesses. In the off-mass shell matrix element the virtual gluon polarization tensor is 
taken in the form of the SHA prescription 0, 0 : 

= ^2^^ = pVa;V|grr = QtQt/\Qt\‘^ (2) 


Nowadays, the signihcance of the kT factorization (semihard) approach becomes more 
and more commonly recognized. Its applications to a variety of photo-, lepto- and hadropro- 

i-El 


duction processes are widely discussed in the literature 


Remarkable agreement is 
and elec- 


12 


found between the data and the theoretical calculations regarding the photo 
troproduction |^, [T^ of D* mesons, forward jets [|T^ [T^, as well as for specihc kinematic 
correlations observed in the associated D*+jets photoproduction 
hadropro duct ion of beauty |^, 0, Xc and J/T ^ at the Tevatron. 


17| at HERA and also the 
The theoretical 
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Figure 1: QCD diagram for the open heavy quark production. 


predictions made in ref. has triggered a dedicated experimental analysis of the J/'h 
polarization (i.e., spin alignement) at HERA conditions. 

Here study of heavy quark production processes in pp—collisions at TEVATRON with 
the SHA is presented. Similar investigations have been done by many authors 1,11,11 
earlier. But in the recent papers HU 


obtained results contradict each other. 


We investigate the sensitivity of the inelastic cross section of heavy quark production 
processes at TEVATRON to different uintegrated gluon distributions. Special attention is 
given to the unintegrated gluon distributions obtained from BFKL evolution equation. The 
outline of our paper is the following. In sect. 2, we give the formulas for the cross sections 
of heavy quark hadroproduction in the SHA of QCD. Then, in sect. 3, we describe the 
unintegrated distributions which we use for our calculations. In sect. 4 we give the matrix 
elements for gluon-gluon fusion QCD subprocesses. In sect. 5, we present the results of our 
calculations. Finally, in sect. 6, we give some conclusions^. 


2 SHA QCD Cross-Section for Heavy Quark Produc¬ 
tion in pp Collisions 

We calculate the total and differential cross sections of heavy quark production pp —>■ QQX 
via the gluon-gluon fusion QCD subprocess (Fig.l) in the framework of the SHA. First of 
all we take into account the transverse momentum of initial gluons (^i, 2 r), their virtualities 
(?i 2 = ~^ 2 t) the alignment of gluon polarization vectors along their transverse mo¬ 
menta given by (2) il, 1,0. Let us dehne Sudakov variables of the process pp —>• QQX 

^The obtained results have been reported at Workshop DIS 2001, Bologna, 27 April - 1 May 2001. 
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(Fig.l): 


Pi = ttiPi + I 3 iP 2 + Pit, P2 = ^2^1 + P2P2 + P2T 

qi = XiPi + qiT, q 2 = X 2 P 2 + q 2 T, (3) 


where 

2 271 ^2 2 2 2 2 

Pi — P2 ~ ^ 1 Ql ~ QiTi ^2 ~ Q 2T1 

Pi, p 2 and gi, q 2 are 4-momenta of the heavy quarks and the gluons respectively, pir, P 2 T 
and qiT, q 2 T are transverse 4-momenta of quarks and gluons. In the center of mass frame 
of colliding particles we can write Pi = (P,0,0,P), P 2 = (P,0,0, —P), where E = 

Pi = P^ = ^ and (P 1 P 2 ) = s/2. Sudakov variables are expressed as follows: 


Mir Mot 

tti = ^^exp(pi), a2 = 

ijs yjs 

'T Ado'T 

Pi = ^^exp(-pi), /32 = ^^exp(-p;). 






(4) 


where M/ 2 T — +Pi 2 T 7 Hi 2 rapidities of heavy quarks, M is heavy quark mass. From 
conservation laws we can easly obtain the following relations: 

qiT + q 2 T = Pit + P 2 T, xi = ai + ^2, X2 = Pi + P2 (5) 


The differential cross section of heavy quark production has the following form: 

MpP -^QQX) = !!^ 4 .( 11 , 4 ? t . 0 ") ^ dqtr X 
X 4It. «") § iqir Mrs‘ ^ QQ), 


( 6 ) 


where <F(a;, gf., Q"^) is an unintegrated gluon distribution in the proton, da{g*g* QQ) is 
the differential cross section of gluon-gluon fusion subprocess. We used the following form: 


da{g*g* 


o<5) = ^ y iMp 


2s 


SHA 


{ 9 * 9 * - QQ) 


dPpi dPp2 


5 '^^\qi + 92 -Pi -P2), 

(7) 


(27r)^2pi (27r)3 2p2 
where \^\‘sha{9*9* QQ) is the off mass shell matrix element and 
5 (gi + g2 - Pi - P2) = 5 (g? Pq^-Pi-pl)^(P it + q2T - Pit - P2t) d{ql + ql - pi - pi). ( 8 ) 


In (7) X) indicates an averaging over polarizations of initial gluons and a sum over polariza¬ 
tions of hnal quarks. We have also 

d^pi d^P2 _ d^PiT d^P2T , ^ xg. 

(27r)3 2p? (27r)3 2p^ 2(27r)3 2(27r)3 


Integrating out the P 27 -, xi and X 2 dependences in ( 6 ) and (9) with accounting of (7) and ( 8 ) 
we obtain the following formula for the differential cross section of the process pp —QQX 
in the framework of the SHA: 


da{pp ^QQX) = 4>(a;i, , qf^, <5^) ^( 3 ^ 2 , q^T, Q^) x 

X E \M\Ij,A9*9* - QQ) dy*i dy*2 dplT dqlr dqh 
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If we take the limit qiT 0, q 2 T 0 and if we average (10) over the transverse directions 
of the vectors qiT and we obtain the formula for the the differential cross section of the 
process pp QQX in the standard parton model; 

da{pp qqX) = Q ) ^ \M\],j^{gg QQ) dyl dy^ dpl^, (11) 

where 1^1 pm (S'S' QQ) is matrix elements of gluon-gluon fusion QCD subprocess in 
the standard parton model (SPM). Here idicates averaging over polarizations of on-shell 
initial gluons and a sum over polarizations of hnal quarks. We average over the transverse 
directions of and using the following expression: 

SdiRr <i1t) Sdrir qlr) - QQ = , ,, 

= xiG(xi, Q^)x 2 G(x 2 , Q^) E |A'/|pm(99 ^ QQ), 


where 


2w 


0 


d(j)l,2 Qi,2T Qi,2T 
27r ^,2T 



(13) 


3 Unintegrated Gluon Distribution Functions 


Various parametrizations of the unintegrated gluon distribution used in calculations are 
discussed below. First, as in the publication Bi, we used the following phenomenological 
parametrization (LRSS-parametrization) [^, §]: 


where 


<1>(X, (fp^) = ^’c 


0.05 


0.05 -I- X 


(1 - xffi{x,qT^), 


fi{x,qT^) 


1, iiqT^<qi{x), 

{ql{x)/qT^)\ \iqT^>ql{x) 


(14) 


(15) 


with ql{x) = ql + exp(3.56yZn(xo/x)), ql = 2 GeV^, A = 56 MeV, Xq = 1/3. The 

value of the parameter qQ{x) can be considered as a typical transverse momentum of the 
partons in the parton cascade which leads to natural infrared cut-off in semihard approach. 
The normalization factor of the structure function <l)(x, gp ^)($o = 0.97 mb) was obtained 
in P, from the beauty production in pp at SppS energyf]. The effective gluon distribution 
xG{x,Q‘^), which was obtained from eq. (1) with (14) and (15), increases at not very 
low x(0.01 < X < 0.15) as xG{x,Q^) ~ x~^, where A Ri 0.5 corresponds to the QCD 
pomeron singularity given by summation of leading logarithmic contributions (agln-)^ 
This increase continues up to x = Xq, where Xq is a solution of the equation go(^o) = Q^■ In 
the region x < Xq, there is saturation of the gluon distribution function: xG(x, Q^) ~ 
(curves 3 in Fig.2). 

^The experimental data relate to relatively low energy for applicability of the SHA. 
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Another parametrization is based on the nnmerical solntion of the BFKL evolution equa¬ 
tions (RS-parametrization), which has the following form PB|: 




.“h K + + 

4v 2 TT-^ 02 as 04 q^ 


ax -\- 


p 


]C, 


05 


11/2 


e-I-ln(l/a:)'^ ^^Os-l-x 


[1 


060:“^ In (gVa8)](l + aiia:)(l - 


where 


Cn = 


1 , < qo{x), 

loM/Y, 'PY > loM- 


(16) 


(17) 


All parameters (see [^) (oi — an, a, (3 and e) were found by minimization of the differ¬ 
ences between left hand and right-hand of the BFKL-type equation for unintegrated gluon 
distribution <F(x, at Ql =4 GeV^. 

Then we use the results of a BFKL-like parameterization^ of the unintegrated gluon 
distribution $(a;, g|^,/i^), according to the prescription given in |^. The proposed method 
lies upon a straightforward perturbative solution of the BFKL equation where the collinear 
gluon density xG{x,fi‘^) from the standard GRV set is used as the boundary condition 
in the integral (1). Technically, the unintegrated gluon density is calculated as a convolution 


of collinear gluon density G{x,^^) with universal weight factors pbj : 


4>(x, ql, ^(r/, q^, - G(-, d^, 

Jx 7j rj 


where 


^ Jo{2^as\n{l/7]) ln(p7g|,)), 

mr 


ln(l/? 7 ) ln(gf.//i2)), 


mx 


2 ^ 2 

«T > 1* . 


(18) 


(19) 


( 20 ) 


where Jq and Jq stand for Bessel functions (of real and imaginary arguments, respectively), 
and dig = ScKs/tt. The parameter dg is connected with the Pomeron trajectory intercept: 
A = Q;541n2 in the LO and A = Q;541n2 — Nal in the NLO approximations, respectively, 
where iV ~ 18 |^. The latter value of A have dramatic consequences for high energy 
phenomenology However some resummation procedures proposed in the last years lead 


to positive value of A(~ 0.2 — 0.3) |^, |^. Therefore in our calculations with (18) we used 
only the solution of LO BFKL equation and considered A as a free parameter varying it from 
0.166 to 0.53. Pomeron intercept parameter A = 0.35 was obtained from the description of 
Pt spectrum of D* meson electroproduction at HERA |^. We used this value of the the 
parameter A in present paper. 

Finally we tried the so called differential unintegrated gluon structure function of a 
proton, which could be obtained by differentiation of the collinear gluon density xG{x, Q^) 0, 
for example one from the standard GRV set according to (1) 


Of course LRSS and RS parametrizations are BFKL - type too. 
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Figure 2: Effective gluon densities xG{x, Q^) obtained from different unintegrated gluon 
distribution. Curve 1 correspond to the GRV collinear gluon density [^], curves 2, 3, 4, 5 
correspond to RS (at Qg = 4 GeV^), LRSS (at Ql = 2 GeV^), BFKL @] (at Qg = 1 


GeV^) parametrizations and differential (at Qg = 1 GeV^) gluon distribution. 
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Figure 3: 


QCD diagramms of gluon-gluon fusion subprocesses for heavy quark production. 




dxG{x, 


( 21 ) 


d In Q2 

In expression (1) the accounting of the contribution of the unintegrated gluon distribution 
to integral at low region (0 < gf. < Ql) was done by the collinear gluon density xG{x, Ql), 
where = (1 -i-4) GeV^ in dependence on different parametrizations of unintegrated gluon 
distribution. 

Fig. 2 shows the x dependence of the effective gluon density xG{x, Q^) at hxed val¬ 
ues of obtained with help of the dehnition (1) for different parametrizations of uninte¬ 
grated gluon distribution discussed above. The curve (1) corresponds to the standard GRV 
parametrization of collinear gluon distribution function. Gurve 2, 3, 4, 5 correspond to effec¬ 
tive gluon densities obtained from unintegrated gluon distribution functions <F(a:, q^) with 
help of the dehnition (1) for RS, LRSS, BFKL parametrizations and the dehnition (21) at 
diherent values of Qq parameter. 

The ehective gluon density xG{x, Q^) obtained for the LRSS parametrization of uninte¬ 
grated gluon distribution $(a;, q^) dihers strongly from other gluon densities normalizated to 
GRV collinear gluon density. It is result of the normalization of the LRSS parametrization 
of <h(a;, q^) obtained from the SppS experinental data for bb production in pp collisions 


4 Matrix Elements for Gluon-Gluon Fusion QCD 
Subprocesses 

The subprocess g*g* QQ is described in the LO pQGD by three diagramms shown in Fig. 
3. In according with the standard Feynman rules for QGD diagramms: 
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Pi - qi + m 


Ma = u{pi) i-igj'") e^{qi) i 


{pi — gi)^ — m? 


-igY)£u{q2)v{p2), 


Mb = u{pi) i-igY) £u{q2) ^ ^ (-^97^) 

[pi - q2) ~ 

Me = u{pi) -g2, qi + ^ 2 ) ^ 7a-^^(^ 2 ), (22) 

(gi + q2r 

where Sfj,{qi) and euiq 2 ) are the polarization vectors of the initial glnons, C^'^^(gi, g 2 , gs) 
is the standard QCD three glnon vertex: 

12 , Iz) = Kiiz - + (93 - 92 )'“ g-" + (91 - izYa^n- (23) 

The Mandelstam variables for glnon-glnon fnsion snbprocess g*g* —*■ QQ are 


s = (gi + g2)^ = (pi +P2)^ 
t = {q2-P2Y = (Pi - gl)^ 
u = (gi -P2f = (pi - g2)^ 


(24) 


and 


s i u — 2 T q^rp T q'^j'- 
Let present the matrix elements in the following forms: 


(25) 


Ma,b,c — (26) 

then averaging over initial glnon polarizations and snmming over hnal qnark polarizations 


we have 


where 


E |M|2 = E \M\\ + CbE \M\l + CcE \M\l+ 
-2CabY. \M\\b + 2CacE \M\\c + 2CbcE 


(27) 


^|M|2 = \LfML[%2)MrMr^ 


(28) 


Here the index i = A, B, C, AB, AC, BC and the color factors are Ca = Cb = 1/12, 
Ca = Cb = 1/12, Cc = 3/16, Cab = —1/96 Cac = Cbc = 3/32 accordingly. The effictive 
gluon polarization tensor L^^^q) is taken from (2). 

The calculation of E \MWBAi9*9* QQ) done analitically by REDUCE system. 
The obtained results coincide with the ones from ref. [^]. 

For the calculations of E \M\pj^{gg QQ) in the framework of the SPM we used the 
following form of the gluon polarization tensor in axial gauge: 


92 ) — -g^'' + 


qU2+9U^2 


- 9 I 


qUi 


{qi.q2) ^2)2’ 

Obtained expression for E \M\pj,^{gg —>■ QQ) coincides with the results of ref. p3 . 


(29) 
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Figure 4: 
in Fig. 2. 


The bb production cross section at ^/s = 1.8 TeV. Curves 1-5 are the same as 


5 Results of Calculations 


The calculations of the heavy quark hadropro duct ion cross section in the SHA have been 
made acceding to eqn. (10) for > Qq GeV^. For qr^ ^ Qo GeV^ we set |^| = 0 in the 
matrix elements of subprocesses, take QQ) instead of \M\ %ha{9*9* QQ) 

and use the eqn. (11) of the standard parton model (SPM). The choice of the critical valne 
of parameter gp ^ = gQ = (1 4) GeV^ is determined by the reqnirement of the small value 

of asiqr'^) in the region gp^ > (1 -h 4) GeV^, where in fact as{qT^) < 0.26. 

The results of our calculations for the cross sections of bb prodnction at TEAVTRON are 
shown in Figs. 4 8. These hgures show the dependence of the bb prodnction cross 

section at TEVATRON energies calculated with mb = 4.75 GeV, Aqcd = 150 MeV and 
\yl\ < 1 ( Fig. 4) [|], and \yl\ < 1, \y^\ < 1, par > 6.5 GeV ( Fig. 5) [|] [|. In Fig. 4, 5 
the cnrves 1 correspond to the calculations in the framework of the SPM using the gluon 
density xG{x,Q‘^) from the standard GRV set The curves 2, 3, 4, 5 are the results of 


calculations in the framework of the SHA using the unintegrated glnon distribntion in RS 
(Qo = 4 GeV^), LRSS {Ql = 2 GeV^), BFKL {Ql = 1 GeV^) parametrizations and the one 
(22) {Ql = 1 GeV^) accordingly. 

One can see that the SHA curves 2, 4 and 5 describe the DO and GDF experimental data 
very well in comparison with the SPM results (curves 1). The LRRS parameterization of the 
nnintegrated glnon distribntion overestimates the experimental bb prodnction cross section 

ii- 

We tried to describe these experimental data with help the LRSS paremetrization htting 
the Ql, mb and Aqcd parameters in region: 1 GeV^ < Qo ^ 4 GeV^, 4.5 GeV < mb < 5.0 
GeV and 100 MeV < Aqcd ^ 250 MeV. We obtained well agreement of the theoretical 

^The cross section does not depend on region 1 GeV< p™” — ^ GeV approximately. 
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Figure 5: The bb production cross section at y/s = 1.8 TeV. Curves 1-5 are the same as 

in Fig. 2. 



Figure 6: The bb production cross secton at y/s = 1.8 TeV. Curves 3 and 3* correspond 
to the SHA calculations with the LRSS unintegrated gluon distribution at nib = 4.75 GeV, 
= 150 MeV, Ql = 2 GeV^ and nib = 5.0 GeV, Aqcd = 100 MeV, Ql = 1 GeV^ 
accordingly. 
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Figure 7: The bb production cross secton at a/s = 1.8 TeV. Curves 3 and 3* are the same 
as in Fig. 6. 



Figure 8: The bb production cross secton at a/s = 1.8 TeV. Curves 2 and 2* correspond to 
the SHA calculations with the RS parametrization of unintegrated gluon distribution using 
the matrix elements (22) and the ones from ref. |^. 
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results (curves 3* in Figs. 6 and 7) with the TEVATRON experimantal data at the following 
values of the parameters: Ql = 4 GeV^, mb = 5.0 GeV and Aqcd = 100 MeV [|. 

In Fig. 8 we show the results obtained using the off mass shell matrix elements of gluon- 
gluon fusion subprocesses in the form (22) and the ones from ref. [Q with using the RS 
parametrization of unintegrated gluon distribution as example. We see that these results 
coincide very well. 


6 Conclusions 


We considered the process of inelastic heavy quark production at TEVATRON in the frame¬ 
work of the semihard QGD approach with emphasis on the BFKL dynamics of gluon distribu¬ 
tions. We investigated the cross section a{pT > p™”) of inelastic 6—quark hadropro duct ion 
as a function of different unintegrated gluon distributions. It is shown that the description of 
the 6—quark inelastic cross section at TEVATRON energies is achieved^ in the cases of the 
RS, BFKL parameterizations and also the parametrization (22) at realistic values of QGD 
parameters (mb = 4.75 GeV and Aqcd = 150 MeV). The LRSS parametrization describes 
the DO and GDF experimental data with another values of parameters ( mb = 5.0 GeV and 
Aqcd = 100 MeV) in comparison with the description of the bb quark photoproduction data 
at HERA in the framework of the SHA |^. 

One of the main goals of our investigations consist of a search in the framework of the 
SHA of the ’’universal” unintegrated gluon distribution. In our opinion the studies of many 
high energy heavy quark production processes in the framework of the SHA BIT^, 1^ |T^ ^ 


have showed that the so called BFKL unintegrated gluon distribution function ^6 
the candidates for a role of universal gluon distribution. 


is one of 
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